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In general relativity, systems of spinning classical particles are implemented into the canonical 
formalism of Arnowitt, Deser, and Misner The implementation is made with the aid of a 
symmetric stress-energy tensor and not a 4-dimensional covariant action functional. The formalism 
is valid to terms linear in the single spin variables and up to and including the next-to-leading order 
approximation in the gravitational spin-interaction part. The field-source terms for the spinning 
particles occurring in the Hamiltonian are obtained from their expressions in Minkowski space with 
canonical variables through 3-dimensional covariant generalizations as well as from a suitable shift 
of projections of the curved spacetime stress-energy tensor originally given within covariant spin 
supplementary conditions. The applied coordinate conditions are the generalized isotropic ones 
introduced by Arnowitt, Deser, and Misner. As applications, the Hamiltonian of two spinning 
compact bodies with next-to-leading order gravitational spin-orbit coupling, recently obtained by 
Damour, Jaranowski, and Schafer 4], is rederived and the derivation of the next-to- leading order 
gravitational spin(l)-spin(2) Hamiltonian, shown for the first time in Q], is presented. 

PACS numbers: 04.25.-g, 04.25. Nx 
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Full implementation into canonical formalisms of general relativity (GR) and applications have so far found only 
(SI ' classical point masses [H, i, i, 0, i, i, 113, jlMi, [ll Q [H , fluids [II [lE [11 [M lip , massive scalar fields [H], [p , 
J> ' and gauge spin-1 fields, including Maxwell [l|, |22|| and Yang-Mills [l^. The canonical implementation of spin-i Dirac 
^^O , fields has been performed in [2J, [2^, [53, l27l. l28l|. Formally showing derivative coupling to the spacetime metric, 
• the Dirac fields resemble to the classical spinning objects (pole-dipole particles) treated in our paper. Problems of 
canonical gravity with derivative-coupled sources are discussed in the comprehensive review by Isenberg and Nester 
[29I I. Another common feature of Dirac fields and classical spinning objects is the occurrence of surface terms in 
' the Minkowski space algebra of the stress-energy tensor, see Appendix |^ The canonical formulation of Dirac fields 
I coupled to gravity is therefore a valuable guide for the considerations in this paper and for future work. 
' Regarding classical pole-dipole particles in GR, see, e.g., [sO, [Ml, [12] , the theory of special relativity (SR) tells us 
that only for specific spin supplementary conditions (SSC), namely, the Newton- Wigner (NW) ones [33], canonical 
variables can be achieved. Related with a SSC is an implicit association of the used coordinates to a specific center 
for the particle: in case of the canonical NW SSC the center is called center-of-spin, in case of the noncanonical 
5^ ■ noncovariant M0ller (in SR) or Corinaldesi-Papapetrou (in GR) SSC [H, HH] center-of-mass or center-of-energy, and 
a I in case of the covariant Fokker-Synge-Pryce (in SR) or Tulczyjew (in GR) SSC, [HIH, [13, , center-of- inertia, see, 
e.g., ^3^. If one is interested in a theory with terms linear in spin only, the Fokker-Synge-Pryce- Tulczyjew SSC are 
identical with the Lanczos (in SR) [13] or Mathisson and Pirani (in GR) SSC, [4l|, |^], for history see, e.g., fisj . 

In this paper, the canonical formalism by Arnowitt, Deser, and Misner (ADM), see [l|, will be applied to put the 
GR-dynamics of pole-dipole particles into canonical form. The starting point will not be a covariant action functional 
but rather the symmetric stress-energy tensor of pole-dipole particles. The developed formalism is valid to terms linear 
in spin and, in post-Newtonian framework, to next-to-leading order approximation in the spin interaction part. The 
formalism is applied to the derivations of the ADM Hamiltonian of two spinning compact bodies with next-to-leading 
order gravitational spin-orbit coupling recently obtained in and to the calculation of the next-to-leading order 
gravitational spin(l)-spin(2) Hamiltonian. The outcome of the latter calculation has been announced in [3] already. 
It is hoped to develop the canonical formalism to higher orders in future. 

The canonical dynamics is only given in a reduced form in this paper, i.e. all gauge degrees of freedom due to general 
coordinate invariance are already fixed. A similar reduced formulation for gravitating Dirac fields is given in [25| . 
Thegauge independent canonical formalism for Dirac fields is achieved for tetrad gravity instead of metric gravity 
in [28[, i.e. the vierbein instead of the metric is the fundamental dynamical variable. An analogous canonical theory 
for classical spinning objects would be very desirable. Of course, other methods are also well suited to incorporate 
spin effects into the post- Newtonian expansion of general relativity [U HE, i46i] . However, a common problem of all 
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formalisms, if one aims at a Hamiltonian formulation, is to get Poisson brackets for the variables, or to find variables 
that allow for standard Poisson brackets. Here the ADM formalism presents itself as valuable because one is always 
close to the exact canonical formulation of point-masses and the connection to the global Poincare algebra has already 
been studied in detail in the literature, see, e.g., [13, H^. The global Poincare algebra seems to be the smartest tool 
to construct or validate Poisson brackets within a post-Newtonian setting. Regarding interaction terms nonlinear in 
spin, the ADM formalism has been proven useful too. Beyond leading order, various new non-linear-in-spin binary 
Hamiltonians have been derived recently, [ioj . For sake of completeness it should be mentioned that in [soj a covariant 
action functional approach to the dynamics of pole-dipole particles in external gravitational fields has been introduced 
in Routhian form using vierbein fields and in [5l| the same dynamics has been treated within the language of forms. 
A Lagrangian approach is presented in (s^ . In neither of the latter cases dynamical canonical gravity has been 
envisaged. 

It is important to point out that we will count post-Newtonian orders, i.e. orders in c~^, only in terms of velocity of 
light c originally present in the Einstein field equations. Then both linear momentum and spin are counted of the order 
zero. The next-to-leading order in the spin interaction part therefore appears at the second post-Newtonian order in 
this paper. This makes no statement about the numerical value of these contributions, which can, of course, be much 
smaller compared to the second post-Newtonian point-mass contributions (depending on the numerical values of the 
spin variables). Some papers already respect in their post-Newtonian expansion that the numerical value of the spin 
variables is assumed to be of the order c~^. Then the next-to-leading order spin-orbit and spin-spin contributions, 
both second post-Newtonian in our way of counting, are referred to as second-and-a-half and third post-Newtonian 
contributions, respectively. 

The paper is organized as follows. In Sec. II, the structure of the ADM formalism is outlined. Emphasis is put 
on the role the stress-energy tensor of the matter source of the Einstein field equations plays in the ADM formalism. 
In Sec. Ill, the matter Hamiltonian and its relation to the covariant 3-space components of the matter stress-energy 
tensor are discussed. The Sec. IV is devoted to the stress-energy tensor of pole-dipole particles in Minkowski space 
in canonical variables. The components of the stress-energy tensor occurring in the curved spacetime Hamiltonian 
are constructed by 3-dimensional covariant generalization. The canonical linear momentum is identified as generator 
of the global Poincare algebra. The action functional for centcr-of-mass and spin motions is given. The Sec. V shows 
how the components of the stress-energy tensor in the Hamiltonian can be directly obtained in curved spacetime. 
In Sec. VI, consistency of the obtained formalism is proved to the approximation of the Einstein field equations 
treated in the paper. The Sec. VII is devoted to applications. The next-to-leading order gravitational spin-orbit 
and spin(l)-spin(2) Hamiltonians are calculated. In Sec. VIII, an independent derivation of the next-to-leading order 
gravitational spin(l)-spin(2) Hamiltonian is given using the lapse and shift functions which are not involved in the 
calculation of the ADM Hamiltonians of Sec. VII. Finally in Sec. IX, the Poincare algebra is shown to hold to the 
order of approximation of the developed formalism. The Appendix A presents the local stress-energy tensor algebra 
for pole-dipole particles in Minkowski space and the Appendix B gives the local stress-energy tensor algebra in curved 
spacetime for nonspinning particles. The Appendix C shows the applied regularization techniques. 

Our units are c = 1 and G = 1, where G is the Newtonian gravitational constant. Greek indices will run over 
0, 1, 2, 3, Latin over 1, 2, 3. For the signature of spacetime we choose +2. The short-cut notation ab {— a^^h^ = a^bf^) 
for the scalar product of two vectors and b^^ will be used. Round brackets denote index symmetrization, i.e., 
a^f^b"^^ = ^{a^b'^ + a'^b^). The spatial part of a 4-vector x is x. 



Crucial to the ADM formalism [l| is the Hamiltonian which generates the full Einstein field equations, both the 
four constraint equations and the 12 first order evolution equations for the 3-metric 7^ and its canonical conjugate 



where N and iV* denote lapse and shift functions, which are merely Lagrange multipliers. The super-Hamiltonian 7Y 
and the supermomentum Tii densities decompose into gravitational field and matter parts as follows. 



II. STRUCTURE OF THE ADM FORMALISM 
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Here 7 is the determinant of the 3-metric 7ij = gij of the spacehke hypersurfaces t = const., whereas the determinant 
of the 4-dim. metric g^j^ wiU be denoted g. The canonical conjugate to 7^ is jg^fTr*-' . R is the Ricci-scalar of the 
spacehke hypersurfaces and ; denotes the 3-dim. covariant derivative. The expressions for lapse and shift functions 
are then N — (—5°")"^^^ and N'' = 7'-'goj. For simplicity, we will assume that the relation between field momentum 
TT*-' and extrinsic curvature Kij is the same as in the vacuum case: 

= - Y'l''')Kki (2.4) 

The energy E is defined by 

where , denotes partial space derivatives and (Psi the 2-dim. spatial volume element at spatial infinity. The surface 
expression E makes the Hamilton variational principle well-defined also for variations which do not have compact 
support. After imposing coordinate conditions and constraint equations, 

E turns into the ADM Hamiltonian i^ADM- Comparing these constraints with the Einstein equations, projected onto 
the spacelike hypersurfaces, results in 

^matter ^ ^T'^'n^U, = N ^T^^" , T^^attcr ^ -^Tj'n, = , (2.7) 

where yJ—gT^^ is the stress-energy tensor density of the matter system. The timelike unit 4- vector = [—^^ 0, 0, 0) 
points orthogonal to the spacelike hypersurfaces. The evolution equations of the field, before imposing constraints 
and coordinate conditions, read 

167r dt S-fij ' 167r dt Stt*^ ' ^ ' ^ 

The coordinate conditions must be preserved under this time evolution. These additional constraints fixate lapse and 
shift functions. 

The ADMTT gauge pj, being the most often used and best adapted coordinate condition for explicit calculations, 
is given by 

7u = -I- hj^ , or 37jjj - 7^^,, = , and tt" = . (2.9) 

Here hj^^ has the properties hj^^ = hjj^^ = (transverse, traceless). After imposing the constraint equations (|2.6p . 
the remaining 4 (reduced) field equations read 

1 dnr^j, _ SHadm 1 dhj^ ^ SHadm -^^^ 



167r dt 5hJ^ ' 167r dt Snl^^ ' 

where Tr^rp denotes the transverse traceless part of tt'-' and the variational derivatives must include a projection onto 
the transverse traceless part. We will call the phase space consisting of hj^"^ , yi^Tri^rj, and canonical matter variables 

the reduced phase space, whereas the nonreduced phase space consists of 7ij, j^tt*-' and canonical matter variables. 

The fundamental problem to be solved are the forms of the super-Hamiltonian and supermomentum densities for 
pole-dipole particles in canonical variables. Our approach will be as follows. We first construct the stress-energy 
tensor in Minkowski space with canonical variables. Then taking into account that, respectively, 7^™^*'°'' and 7^™''"°'' 
are scalar and covariant vector densities with respect to 3-dim. coordinate transformations, we put these expressions 
into 3-dim. covariant forms (this procedure had been suggested already by Boulware and Deser [IJl). Afterwards 
we show that the same result can be obtained by Lie-shifting certain components of the stress-energy tensor with 
rotational-free parallel transport of the linear momentum fields. 



III. CONSISTENCY CONDITIONS 



The Hamilton variational principle must generate the Einstein equations. This trivial fact leads to several consis- 
tency conditions for the matter part of the Hamiltonian, 



^matter ^ / (^^^.(^T^mattor _ ^»-^mattor) _ (^3 j^-j 
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Lapse and shift are Lagrange multipliers, so Ti™'^**®'' and 7^™'^**'^'' rnust be independent of them. Equation (|2.7|) 
then already ensures that the constraint Eqs. (|2.6p are correct. The evolution Eqs. (|2.8p coincide with the Einstein 
equations if and only if 

r rrinatter 

^ = »- (") 

r rrmatter 

= 2^^^-- • (3.3) 

Violation of the first condition would give an incorrect evolution equation for . This is critical, because the geometric 
meaning of this equation is the definition of the extrinsic curvature Kij = —NT'^^ , additional terms here imply leaving 



Riemannian geometry. This might be fixed by adjusting Eq. (|2.4|) . see 55|. The second condition ensures that the 
evolution equation for tt*-' fits with the Einstein field equations. 
The first condition, Eq. (|3.2p . is equivalent to the local equations 

The second condition, Eq. I|3.3p . then implies that also Ty is independent of tt*-'. If and only if Tij does also not 
depend on lapse and shift, then the local version of the second condition reads 



In Appendix |B] it will be shown that Eq. (|3.5p is equivalent to the simple constraint algebra (jBlHB3|l . The conditions 
given in Eq. (|3.5p are very restrictive, as they imply that 7^™^**°"^ cannot depend on derivatives of , and 7^™^**°"^ 
cannot depend on 7*-' at all. Together with ()2.4p . this defines a kind of simple coupling of matter to gravity. Gravitating 
classical spinning objects and Dirac fields arc not of this kind. However, our canonical formulation of spinning objects 
will exactly fulfill (|2.4p and (|3.4p . and also at least approximately (|3.3p . see Sec. IVII 

None of the preceding consistency conditions validates the canonical matter variables directly, in our case position, 
linear momentum and spin. In a theory that is of the simple kind mentioned above, this can be done via a local 
algebra for 7^™'***'='' and 7^™'^**'"' on the nonreduced phase space, Eqs. (jBllllBT3| . We will instead consider the global 
Poincare algebra, which is a consequence of the asymptotic flatness and is represented by Poisson-brackets of the 
corresponding conserved quantities. So besides the ADM energy (|2.5p also total linear momentum Pi, total angular 
momentum Jj = ^eijkJjk and the boost vector are conserved and given by surface integrals at spatial infinity. 
The boosts have an explicit dependence on the time t and can be decomposed as = G' — tPi, where X' = /E 
is the coordinate of the center-of-mass. G" will be called center-of-mass vector in the following. The corresponding 
surface integrals read, with spatial coordinates denoted x^: 

P^ = ^d^SfcTT^'^ , = -i- (j^ dhkix'TT^'' - TT^'') , (3.6) 

After imposing constraints and coordinate conditions, these quantities have well-defined Poisson-brackets on the 
reduced phase-space [13], and the Poincare algebra can be verified. At the second post-Newtonian level for spin, and 
also in the spatial conformally flat case 7^ = 'ip^Sij, we have, by virtue of the momentum constraints Tii = and the 
ADMTT gauge, the following simple expressions: 

P,^ J d^xnT''''°' , ^ J d^x{x'nf ''''''' -x^HT'"''''')- (3.8) 

In the ADMTT gauge it also holds: 

E = Hadm = -^ I d^xAi;, I d^xx'Ai;. (3.9) 

After solving the Hamilton constraint Ti = 0, tp can be expressed in terms of canonical variables of the reduced phase 
space, and the Poincare algebra can be verifled. 

A flnal remark concerns the canonical spin variables. Imposing the standard Poisson-bracket algebra of angular 
momentum for the spin variables, it is clear that the squared euclidean length of the spin, being a Casimir operator, 
will commute with all other canonical variables. Therefore this length is constant in time, as it will also commute 
with the Hamiltonian. 
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IV. POLE-DIPOLE PARTICLE STRESS-ENERGY TENSOR IN CANONICAL VARIABLES 



Calculating 7^™''**'"' and T^f via (|2.7|) in the Minkowskian case and then going over to their 3-dim. covariant 
generalizations has the advantage that TY'"***'"'' and 7^™''*'°'' will definitely not depend on lapse, shift, and tt*-' or 
K^^ . This is a serious problem when working in curved spacetime. Then our matter variables (in particular, spin 
and momentum of the particles, but not their position) have to be redefined to suit the consistency conditions of 
the previous section. It was also observed in 54] that the correct general relativistic source terms in the constraint 
equations for low spin (< 1) fields, including electrodynamics, can be achieved by expressing their flat-space action 
in a 3-dim. covariant form, and redefining canonical variables in a way that leaves them unchanged in the flat case. 
This is similar to our approach. In the next Section we will show that a curved spacetime approach is also possible, 
yielding the same result as in the present Section. 

Because of its importance for later transition to curved spacetime with canonical variables, the stress-energy tensor 
density for an electric charge-free pole-dipole particle in curved spacetime takes the form, to linear order in spin, see, 

e.g., [Hii, 

T'"' ^ j dT [mu^u'^(5(4) - (S'"^''u'''(5(4))||a] (4.1) 

= p^'v''5 ~ (5"(^w'^)(5)^„ - S°'^^'T''^pV^ 5 , (4.2) 

applying the Tulczyjew SSC or, equivalently, the Mathisson-Pirani SSC 

S'^'^u^^O. (4.3) 

Here — jvP and = mu^, particularly pi = mui, with mass m and g^^u'^u'^ = —1. The Christoffel symbols 
are denoted F^^, as usual, and the 4-dim. covariant derivative by ||. The 4-dim. spin tensor S^'^ has the property 
S^'^ — —S'^^. r is a proper time parameter running from — cx3 to + cxd with = where is the 4-dim. 
position variable of the particle. The coordinate time velocity of the particle, , is identical with The Dirac 

delta functions, (5(4) = 5{x — z) and 5 = (5(x — z), are normalized such that / d'*a:(5(4) = / (Px5 = 1 holds. 
Furthermore, again to leading order in spin, it holds 

— = 0, (4.4) 

where D denotes the 4-dim. covariant differential. Obviously, 

S^"'S^„, = 2s^ ^ const. (4.5) 

is valid. 

The transition to Minkowski space results in the stress-energy tensor 

T^"' ^ J dr mu'^u'^(5(4) - (S'"^'' it") (5(4)), 

= p^f''(5- (S'"(^w'')(5),„. (4.6) 
Now we proceed to the Newton- Wigner SSC in making the following shift of the particle coordinates 



m ~ np 



(4.7) 



where np = n/^ip'^ — — -y/ -t- j^^PiPj , as well as introducing the spin tensor S"^" by the relation, see, e.g., [39|, 

gf.^ = ^M"' + pt^nxS^^/m - p^'nxS^'^/m , (4.8) 

which results in 

{n,+p,/m)S^'' = 0. (4.9) 
This turns the stress-energy tensor into the form (from now on S = (5(x — z)) 



f^"'{x, z) = T^"'{x, z) = p^'v''S - {S"'^^'v''^ S) ^. 



(4.10) 
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because z^^ = (dot means time derivative) to linear order in spin. The new spin tensor S^'^ has the important 
property that 

S>"'S^^ = S^jS.j = const. (4.11) 

is vahd. 

The components of the stress-energy tensor, relevant for the ADM formalism, read 

y^f'^'"n,,n^ = -npS - ( dijSki — — — SjkS ] , (4-12) 

-VlTiTLi, = + i (^mkSik - {SmkSip + 6„ip5ik)SqlSqp ^^^'^ ^^-^ ^ • (4-13) 

These components of the stress-energy tensor fulfill the Poisson-bracket algebra a stress-tensor has to fulfill in 
Minkowski space, see [i^,!!^]- Details are given in Appendix [K\ 

The 3-dim. covariant generalizations of these expressions read (; denotes the 3-dim. covariant derivative) 

^matter ^ ^ _j^pS _ /^yj^fc' SjkS] = ^f°° , (4.14) 



m — np 



.^matter _ _ 

Correspondingly, 



-Vlftn^ -P^S+l[ [r'S.'^ - ir'Sf + ^™pjf )y'^,,_^^^) ^ N^f^ . (4.15) 



"f^'^-f^^SijSki = 2s^ = const. (4.16) 



has to hold. The new canonical spin variables 5'(i)(j) (the round brackets make allusion to implicit dreibein components) 
are defined such that 

yS^'^.j^M = = (4.17) 

is valid. This can be achieved by constructing Cij as the symmetric matrix square root of symmetric jij {^ij = ^ji), 

eueij = 7y , Cij = eji . (4-18) 

Then it holds 

Ski = ekteijS(i)(j) . (4.19) 

The condition Cij = Cji had also been imposed on the spatial part of the vierbein field in [25| in order to achieve a 
canonical formalism for the spin-^ field. 
If the 3-metric is represented in the form 

Itj = 5.,j + h,j , 1% I « 1 , (4.20) 
the solution for reads (with some abuse of notation) 



and the variation of e.y is given by 



Sy = \/ 5ij + hij = 6ij + -hij - -hikhkj + ■ ■■ (4-21) 

Z o 



^e^i = ^(57y - ^{hkjS-fik + hikSjkj) + ■ ■ ■ ■ (4.22) 



It may be pointed out that the simple variational relation for dreibein fields e(j)i, where 7^ = ei^k)i^{k)j: of form 
2i5e(j)i = e[j)kl^^5^ii is not valid for our symmetric matrix square root in general; exceptions are isotropic metrics. 
Recalling Eq. (|3.8p . the new canonical momentum Pi is defined in the way that the following structure holds. 



-^matter = + 1 



(4.23) 
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where 

Sjk. (4.24) 



np[m — np) 



^matter ^ _ 1^^. ^ „ ( ^fl^yJ ^kl g^^g^ ^ (4 25) 



Hereof, we get 

^matter ^ _ 

2 ' Vm-nF 
where the quantity can be related to the flat ^/jTij via Eq. (|4.10p : 

nF ■' ■' nF (nF)^[m — nF) 

The crucial question now is for the canonical variables. For both the second post-Newtonian order approximation 
for spin and the spatial conformally flat case in general we get for linear and angular momentum 



(FxTiT = Pi, (4-27) 

Jy = J d^x {x'nf - rr^T^r"") = ^'Pj - PP, + . (4.28) 

It is important that these expressions were achieved in the ADMTT gauge and with — Cji . Under these conditions 
both generators of the global Poincare group flt with the standard Poisson-brackets, 

{z\t),Pj{t)} = S,j , {5(,)(t),5(,)(t)} = e,jkS(k){t) , zero otherwise, (4.29) 

where <S'(i)(j) = ^ijkS^k) {S(i)Si^i^ = s^) with the completely antisymmetric Levi-Civita tensor eij^ (£123 = !)■ In the 
following we will also use the notations S for P for Fi, and Z for z*. The commutation relations of the field 
variables still read 

{hj^{^, i), 7r^V(x', t)} = 167r5™J(x - x') , zero otherwise , (4.30) 

where 

- {Su - A-^dkdi){5,, - A~^did,)] 

with the inverse Laplacian A^^ and the partial space-coordinate derivatives di. Herewith we have completed the 
calculation of the source terms applicable to the ADM formalism. Crucial for our approach is the property of our 
spin variables S to have conserved euclidean length. Further discussion of the consistency of our formalism is given 
in Sec. Ell 

The ADM Hamiltonian, written for a many-particle system (numbering a — 1,2,...) depends on the following 
variables, 

-ffADM = HadM K, -Pai, Sa{i) , hj^^ , TT^t] (4.32) 

and the corresponding action W reads (dot means time derivative) 

^ = J^^^i^ ^"^^ + ^ ^^'^"^'^ ^ i / '^'^'^TT^r - Hadm [zl Pau Si^\hj;^, 4^] ^ , (4.33) 

where n'^^ = \eijkK(i){3)K{i){k) with Ka{i){k)K(j)(k) = K(k)(i)K(k)(j) = Hereof, by variation of W with respect 

to Fai, zl, si''' = ie'^'''S'o(j)(fc), Ao(i)(j) in the forms SFai, 6zl, SSi'\ (SOi*^ = ^eijkAa(i){j)SAa{i){k), the equations of 
motion follow: 

''^^'>- 5Fam ' 5zl{t) ' ^^-""^^ 
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"i'^(^) - ' - ^^^^'^S^it) . (4.35) 

The field evolution is obviously given by Eq. (|2.10p . 
Finally, the transition to a Routhian reads 

R[zl Pa,, Sa(k),hj;^, /i^T] ^ Habm[zI Pa^ , S a^k) , hj^^ , nf^'^] -^jSx n^'^hl'^ , (4.36) 

with the field equations 

5jR{t')dt' _^ 



5hJ^(x\t) 



and the equations of motion 



(4.37) 



The Routhian is very suitable for the derivation of an autonomous, conservative Hamiltonian for the matter, where 
the solution hj^ of the field equations is replaced by the matter variables, see Is^. 



V. SPACETIME APPROACH TO THE STRESS-ENERGY TENSOR IN CANONICAL VARIABLES 

The 3-dim. derivation in the previous section of the needed stress-energy components does not show up which 4-dim. 
object in curved spacetime is behind the performed construction. This will be clarified in this Section. Starting from 
our original curved spacetime stress-energy tensor density with covariant SSC, Eq. (|4.ip . we add up the following 
Lie-shift to it. 



ihiftcd 



dr 



dr 



where = —S'^^n^/{1 — nu) and 



(5.1) 



(5.2) 
(5.3) 



as generalizations of Eqs. (|4.7p . (|4.8p and (14. 9p to curved spacetime. Note that now introduces the lapse function 
into these expressions. Equation (j5.ip was found to be the stress-energy tensor of a spinning particle with mass dipole 
moment in [58| . i.e., its position variable is the Newton- Wigner one in the Minkowski limit. Unfortunately, an 
explicit calculation shows that the components -/V(\/—5r°°) shifted and gti, {^/—gT^'^) shifted still depend on lapse and 
shift, which is not compatible with the ADM formalism. The solution to this problem is inspired by the observation 
that multiplication with and gii, does not commute with taking the Lie-derivative. 

Therefore, we first calculate the projections of the stress-energy tensor density with covariant SSC given by Eq. 



(|4.ip . i.e., y/jT'^^Uf^rLi, and —y/^T^Ui,. These quantities, after a long calculation, turn out to be independent of lapse 
and shift. Adding up their Lie-shifted expressions (notice mP = 0) and also using the definitions (|5.2p and 



Pi 



ik ; 



(5.4) 



which fortunately eliminates Kij and therewith tt'-' , we end up with the expressions 



(V7T^"^n^n,)shifted - -np6 - [ fi^^'—^L^J^^s \ = N^^f"^" 



m — np 



{'-^/lT^n^)shi^tcd = Pi5 + 



1 



PiPk 



2 

-t- Sx^iPu + pi^, - pu)6 = Ny/jTl 



np{m — np) 







(5.5) 
(5.6) 
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where Sx'' = —m^/m. Full agreement is obtained with our previous results (|4.14p and (|4.15|) if the linear momentum 
Pi (as function of space and time coordinates) gets parallel shifted along 5x^ and shows no rotation. Then pi and Pi 
play identical roles and may be identified and thus, T'^" and T'^" too. In order to fulfill the global Poincare algebra, 
we must indeed drop this term proportional to 5x^ . Including it into the definition of our canonical momentum (j4.24p 
is not possible, see Sec. IIXI 



VI. CONSISTENCY CONSIDERATIONS 



Our action (I4.33P has the important properties that it exactly coincides with the expected spin dynamics in the 
Minkowski case, that it reduces to the usual point- mass dynamics for vanishing spins, and that our spin variables 
have constant Euclidean length like in the covariant equations of motion approach for spin, see [2]. Our action, 
formally valid up to arbitrary order, thus defines a spin dynamics that should at least be a good approximation to the 
dynamics described by the covariant stress-energy tensor (|4.ip . We will argue in the following that up to the second 
post-Newtonian order, i.e., the next-to-leading spin-orbit and spin(l)-spin(2) order, our dynamics is indeed the same 
as of the covariant stress-energy tensor treated as source in the Einstein field equations, see, e.g., [4^ . 

First we define 



n — nP 



1 



Then and are simply given by H'"^""' 

(13.51) we now have: 



P P 

nP 



(5(x 



2 

= ~nP5 ^ 
1 ^^^(x) 



nP{m — nP) 
\ and 



Pri 



(6.1) 



Instead of 



/ .rn 



(57^™^""(x) 



2(57*J(x') 

fri(x) 



(x) 



(57*J(x') 



1 



4(x)^(x-x') 



(57^J(x') 



,1 



(57*J'(x') 

At the leading order the total divergences in (|6.2p and (|6.3p do not contribute to (|3.3p 

r rrmattcr i 

:-iVy^4+0(G). 



(57'J 



(6.2) 
(6.3) 

(6.4) 



Note that ^/^Tij is here a Minkowski expression, where our variables are definitely the correct canonical ones. This 
ensures that the evolution equations of hj!^ and Tr^rp are correct at the leading order, which is sufficient for a second 
post-Newtonian Hamiltonian for spin, see also Eqs. (|7.10p and (|7.11[) . 

The structure of (|4.27p and (|4.28p is very promising, as it already implies the fulfillment of a major part of the 
Poincare algebra. This is a very strong argument for our spin variables to be canonical up to the second post- 
Newtonian order for spin and also in the spatial conformally flat case, or, from a different point of view, for our spin 
dynamics (|4.33p to be physical. This argument applies to the reduced phase space in the ADMTT gauge, also recall 
Cij = Cji, where ()4.27p and (|4.28p were derived. The problems encountered for a gauge independent formulation are 
briefly presented in Appendix |b1 

In addition, the next-to- leading order gravitational spin-orbit coupling we will obtain in Sec. I VIII is the same as in 
0. The latter was based on a completely different approach using only the equations of motion for spin; the stress- 
energy tensor for spin was not needed. Also the remaining generator of the Poincare group was determined there 
and the Poincare invariance was shown for the two-body case. In Sec. IIXI we will extend the proof of the Poincare 
invariance to the spin(l)-spin(2) interaction case. 

Finally we present a nice property of our spin variable, both in the second post-Newtonian approximation for spin 
and the spatial conformally flat case. In both cases we can set S'^^/ = Saii"i^^ = Saij^^"^ = Slj = Sa{i){j)- This is 
obvious in the spatial conformally flat case. The neglected hj^'^ contributions are merely total divergences at the 
second post-Newtonian order in the Hamilton constraint, which do not contribute to the corresponding Hamiltonian. 



VII. APPLICATIONS 



In this section we will derive within our formalism the ADM Hamiltonian of two spinning compact bodies with next- 
to-leading order gravitational spin-orbit coupling, recently obtained in and with next-to- leading order gravitational 
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spiii(l)-spin(2) coupling. Some calculations in this and the following Sections were confirmed with the help of xTensor 
[S^ . a free package for Mathematica [60] . 

First we have to solve the constraints iteratively within the post-Newtonian perturbation expansion, which can be 
seen as a formal expansion in c~^. In the source terms of the constraint equations, the action of the mass m has to 
be counted as m ~ 0{Gc^^), and similarly P ~ 0{Gc~^) and S ^ 0{Gc^^). In the following a subscript in round 
brackets denotes the formal order in c~^. We further set ^p = 1 + 0/8 and tt*-' = tt*-' +tt!^j. in our coordinate conditions 
(|2.9p . tt'^ can be written in terms of the vectors tt' = A^^tt*"* j = A^^tt*-' j and tt' = {Sij — ^didjA^^)TT^ as 



(7.1) 
(7.2) 



The Hamilton constraint for an arbitrary source, 



1 



1 2 



to the order needed for a second post-Newtonian Hamiltonian for spin, then reads 



(7.3) 



1 

16^ 
1 



1 

1^ 



a_/mattcr 
rt(2) 



1 

16^ 



'T/mattcr /o/inattcr 
^(6) "oV«(4) 



a_/inattcr /a_/mattcr 



A0(4) = Hflf" 

j -^matter. 



1 



■•(2) 



(4)y 



a_/niattcr j, 
^(2) <P(2) 7 

I -T/mattcr ±2 , ^ 

^64^(2) '^(2) +16^ 



a_/rnattcr j i o/mattcr i \ i /a_/niattcr j2 i oo-zmattcr ^ J. ^ 
(2) -t- At(4) 0(4) -f At(2) 0(6) j + ^(^(4) 0(2) + ^^(2) 0(2)0(4)) 



1 



matter i 3 



1 



g0(2)(^(3))^ + 2^^(3)^(5) - -^<t>{2),^(l>{2),3hf^)^J + ^ (/l(4]'^J,/c)^ 



(td), 



(7.4) 

(7.5) 
(7.6) 



where (td) denotes total divergence, and the momentum constraint. 



jk 



(7.7) 



can be expanded as 



_ o/matter 

167r''(3),.- 2^(3). ' 

I^^5)., - -2^(5). -5^('^(2)^3))..- 



(7.8) 

(7.9) 



The solution to the partial differential equation tt'"' ^ = Att* = A' for tt'-' is given by (|7.2p and tt* = A ^A*. The 
ADM Hamiltonian can now be calculated via Eq. p.9p . 
In the near-zone /i7!,T results from: 



327r- 



""•(4)fc; ^ 



(2),fc<f(2),i 



-167rT(4)M - -0(2),fe0(2)J 



(7.10) 



(7.11) 



The first of these equations is a consequence of the evolution equations (|2.10p , the second is a direct consequence of 
the Einstein equations. Both lead to the same result, if the consistency condition (|3.3p is valid at the leading order. 
At this order Tr^rj, vanishes in the near-zone, the transition to the Routhian (|4.36p is therefore trivial. 

Now we introduce new indices a and 6 that number the spinning particles. Expanding (|4.25p for a many-particle 
system yields 



(2) —2_^^a0a, n 



matter 
(4) 



E 



P2 1 



2ma 2ma 



ai'Ja{i){j) 



5a,] 



(7.12) 
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n-i 
n 



(6) 



'T_/matter 
«(8) 



E 

a 

E 



(P2)2 P2 1 P2 1 



(p2\3 ('P2^2 rp2 p2 1 

16m5 '^(2) a -I V^(2) a zLm 



8m3 



64mQ 



4rna 



(7.13) 

p2 



(td). 



and (|4.23p reads 



(7.14) 



(3)z 



E 



-^matter ^ ^ __ (p^^p^^. ^^^^.^ 4. PajPakSa(j){t)Sa,k) 



(7.15) 
(7.16) 



The leading order of (|4.26p is 

-^(4)u = E 9 C^PaiPajSa + PaiSa(j)(k)^a,k + PajSa{i)(k)^a,k) ■ 



(7.17) 



The equivalence of (fTTOl) and ([7111) can now exphcitly be checked. The source terms of </'(4), ^^3^ and ^^^ij arise 
from the point-mass source-terms by a substitution Pai — > Pai + As this substitution commutes with 

and S]j^^'', we can just apply this substitution to the point-mass solutions of i?!)(4), Tr^g^ and hj^-j, which are, e.g., in 
[11. The results are, with Tq — |x — ^a|; 



^4) 



= 2E 



~ij spin 
^(3) 



, TTspin 



PaiSa(t){j) ( 1 



E 



Sa{k){i 



Sa{k)U) 



E 



PamSa{k](l) 



{^5k(i5j)mdi - 2Sij6kmdi) h (Skmdidjdi - 26k(idjjdmdi)ra 



(7.18) 
(7.19) 
(7.20) 



In order to get this expression for hj^-j, it is actually easier to solve (|7.10p directly, utilizing the formula 87rA ^(5a = 
—Ta, than to use the substitution. The unknown fimctions 0(g) and tt^^^ are not needed for the second post-Newtonian 
Hamiltonian i?2PN = ~tItF / '^'^^^'^(8)1 they disappear after some partial integrations. 0(g) can be eliminated by 



^3 ^ /T/inattcr i 
d X At(2) 0(6) 



1 

16^ 



d'^a:;(A0(2))0(6) = j d''a;0(2)(A0(6)) 



(7.21) 



and then using the constraint (17. 5p for 0(g). Using (|7.2p for Tr^g^ and also ()7.9p and ()7.16p . we get, 

(2)^(3) + E ^^aiP'afe'5'a( 



(3) (5) 



(7.22) 



The /i(4jjj part of the Hamiltonian can also be simplified. We define ^(4)ij such that Ah^-^y^^ — o^^^ /i(4)fc;: 



(7.23) 



The hj^-j contribution to the Hamiltonian then is 



1 

16^ 



^ A f,TT , 1 fiTT \2 
-/l(4)y/l(4).^. + -(,ft(4)y^fej 



1 

16^ 



(7.24) 
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Here we used the fact that Sj^^''^ is a Hermitian operator, (hj^^j j,)^ 
5jj^^^hJ^f,i. The spin part of this can further be written as 



(td), and of course 



1 f j3 ^ ^ point-mass 7 TT spin . ^ ^spin 7 TTspii 



(7.25) 



Note the factor i instead of ;| in the spin-orbit part. This transformation of the hj^^^ contribution is very convenient, 
because the spin part of hj^-^ is much simpler than its point-mass part, which does not contribute any more in Eq. 

The integral H2PN = J d^x A(j)(^g-j can now be computed. The regularization is, at the second post-Newtonian 
order, done by Hadamard's partie finie method and by analytic regularization, see, e.g., [61I [g^. In Appendix [Q 
the formulas needed to regularize the integrals occurring in this calculation are assembled. 



A. Results for Hso° and Hss° 



Now we are ready to present the results. Our Hamiltonian for two spinning compact bodies has a next-to-leading 
order spin-orbit part H^q'^ and a next-to-leading order spin(l)-spin(2) part H^g'-' given by: 



H< 



NLO 
SO 



((Pi X Si) -1112) 



((P2 xSi).ni2) 



((Pi X Si).P2 



((Pi X Si) -1112) 



' 12 



5m2P? , 3(Pi-P2) 



3P? 



3(Pi-ni2)(P2-ni2) , 3(P2-ni2) 



8mf Am\ 4mim2 

' (Pl-P2) ^ 3(PiTli2)(P2Tli2) 

miTO2 mini2 

2(P2 ■ ni 2) 3(Pi-ni2) 

r71lTO2 

llm2 



Ami 



2mim2 



(7.26) 



5m2 
mi 



Anrf 

((P2xSi).ni2) 



' 12 



6mi 



15r7i2 



+ (1^2) 



^SS 



[|((Pi X Si) • ni2)((P2 X S2) • ni2) + 6((P2 x Si) • ni2)((Pi x S2) • ni2) 



2mi'm2ri2 ^ 

- 15(Si • ni2)(S2 • ni2)(Pi • ni2)(P2 • 1112) - 3(Si • ni2)(S2 • ni2)(Pi • P2) 

+ 3(Si • P2)(S2 • ni2)(Pi • ni2) + 3(S2 • Pi)(Si • ni2)(P2 • ni2) + 3(Si • Pi)(S2 • ni2)(P2 • ni2 
+ 3(S2 • P2)(Si • ni2)(Pi • ni2) - i(Si • P2)(S2 • Pi) + (Si • Pi)(S2 • P2) 

- 3(Si • S2)(Pi • ni2)(P2 • ni2) + i(Si • S2)(Pi • P2)] 



(7.27) 



-h((Pi xSi)-ni2)((Pi xS2)-ni2) 



2?71?fl2 

:^4^[-((P2 X S2) • ni2)((P2 X Si) • ni2) 

Zm2Ti2 

6(rni -I- 7712) 



(Si • S2)(Pi • ni2)2 - (Si • ni2)(S2 • Pi)(Pi • 1112)] 
(Si • S2)(P2 • ni2)2 - (S2 • ni2)(Si • P2)(P2 • ni2)] 



[(Si.S2)-2(Si.ni2)(S2-ni2)] 
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Here ri2 = |Zi — Z2I is the euclidean distance between the two particles and ni2 denotes the unit vector ri2ni2 = 
Zi — Z2. (1 <-> 2) stands for repeating the preceding terms with particle one and two exchanged. H^q'^ is identical 
to the result in [2]. The result for H^^^, already announced in Q, differs from the correspondingspin(l)-spin(2) 
potential, V^jj^p^, in [6^. A canonical transformation connecting both results could not be found [J^. In a recent 
preprint , prompted by the preprint version of (3j] , a missing contribution in Eq. (4) of |,63|] has been identified, see 
HI], [Eq. (2)], using information from (gH], [Eq. (18)]. 

The term —j^ J d'^a; 0(2) (vr^gj)^, that contributes to the Hamiltonian via (|7.22p . is the only one where terms 

proportional to Si and S2 survived the regularization procedure. These terms must be dropped, because we already 
neglected them in the stress-energy tensor. 

Of course we are also able to calculate the leading order spin-orbit and spin(l)-spin(2) Hamiltonians via Hipn ~ 
— (167r)~^ J (Px A(j)(^Qy which gives the well-known results: 



s, 



2ma 



(7.28) 
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^sT = ^ E E - ^<^^J (7-29) 

Here r^ft = |Z(j — Z^l and rabn\f^ = — z^. These formulas are even valid for arbitrary many particles. 

VIII. DIFFERENT DERIVATION OF /f|'s^° 

In order to confirm our result for Hgg'~', we use the method from 3] to rederive H^^'-'. Our ansatz for -ffgg^^ linear 
in Si and S2 is now: 

-"SS - ^'(4)y'3l '^2 ^^'(4) -^l— i'(4) • '=>2 • l»-tj 



Note that the equal signs are correct here, because 



--vspin(2) 



^(4)H'5'2''^ (8.2) 



already includes the full dependence of the Hamiltonian on S2. The formula for r2(4) given in 2] can be used without 
further changes, but now the spin-dependent parts of the quantities have to be inserted. The evolution equations, 
correctly given by ()2.8|) if (j3.2|) and ()3.3|) are fulfilled, read: 



7y, - 2N-f-^^\7r,j - hy-Tfci^''') + N,,, + N,,, , (8.3) 

Here R*-* is the 3-dim. Ricci-tensor. Now we determine lapse and shift by demanding that our coordinate conditions 
(|2.9p are preserved under this time evolution. In particular we insert (|8.3p into 2>^ijfij — Ijj.oi — 0, and we take the 
(5jj -trace of (18.41) . The post-Newtonian expansion of the resulting expressions, with further simplifications using the 
constraints, leads to: 

-\<t>i2), (8.5) 

matter _a_/m.atter / \ ^ f J. A \ c;\ 
(4) -^^(2) 0(2) + Y^(0(2)0(2),»),i , (8.6) 

(8.7) 

'^(2)^(3) + ^(3)(i0(2),») . ~ ^[^(3)i<^^>(2)J],^ ■ (8.8) 

Note that also is needed for N. The solution of AiV^ + ^Njji = is given by Ni = ((5,y - ia^^j A"^) A^M^ . 
Again we can get A'^(-4') and A'^(3)i by the substitution Pai + 5'5'(i)(j)5j from their point-mass solutions. This 

gives: 







= 1 , ^(2) = - 




AA^(4) 


= 47rT(4),, -1- 47rH 




1 


= 16717^^3^-, 









J a ^ J 



A'^l^s)^ is more complicated, but for f2(4) we only need 



TvT-spin X ^ 



-Pas-Pam'5'a(m)(i) , r. /^^ 

"2 -2 — I -f maba(j)(m) 

'"■a \'a/ kl \'a/ ,k 



(8.10) 



where Sab = ra + rb + Tab and df and d\ are partial derivatives with respect to Za and Zft, and we used the formula 
A In Safe = (^a''b)^^. Finally, we get from the leading order spin-orbit Hamiltonian (j7.28p : 

^{3)a =Ua,Hso} = ~ + ^^b^j) j ■ (8.11) 



b^a 



ab 
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Now $^(4^2 can be calculated by applying partie finie regularization, e.g., 



where a = 1 and 6 = 2, or a = 2 and b — 1, and Reg^(/(x)) — /i.og(Za), see Appendix [Cl Although this term is 
not symmetric under exchange of both particles, the final result (|8.ip recovers this symmetry, and indeed turns out 
to be the same as (|7.27p . It should be stressed that this approach is indeed independent from the one of the last 
section, in particular, lapse and shift functions had to be determined, also using T^j, and $1(4) was determined using 
the equations of motion of a spinning body in 



, (8.12) 



IX. APPROXIMATE POINCARE ALGEBRA 



At last, the Poincare invariance at the next-to-leading spin(l)-spin(2) order was not yet verified. First we calculate 
Gf^'^ and G^g^^ with the help of jHS]), i.e., G2PN = --^^ J (fxKAcjj^^e)- Using the 3-particle integrals from Ref. ^ 



nrc 



An 




= -47r 









4^LA-i(^^b)]..z.--|^[io.L 



J x=Zc 

and treating the origin as a particle coordinate, results in 

p2 



5r 



ab 



^"l^bc] '''be ) 



(9.1) 
(9.2) 



G 



NLO 
SO 



a " a b^a 



-5{Pa X S,) 



,^NLO 



51:1: 



^(Pb X S„) - i(n,fc X S,)(Pfc • n,b) - ((Pb x S„ 



((Pa X Sa) • liab) 
nab) 



5Za 



Tab 

Za + Zf, 

''ah 



a b^a 



(Sb • nab)% (3(Sa • nab)(Sb • n„b) - (S„ • Sb)) ^ 



(9.3) 



(9.4) 



We get the same result for GgQ*-* as in 0], if we consider our expression for two particles. Now the Poincare algebra 

rNL" ' ' 

^ss 



for two bodies, including H^^'^, can be verified in the same way as in [2], also see [66|. The algebra is indeed fulfilled. 



If we would have kept the term proportional to (5a;' in Sec. |Vl then we must include it into the definition of our 

p2 

canonical momentum (|4.24p . This gives only a change in 7^™^"°'', in particular, the term — g^PaiSa{i)(j)(t>{2),jSa in 
Eq. (Trul) disappears. Then we must add ((^i ^^^i)'"!^) +(1^2) to the Hamiltonian in Eq. FlM . but the 

center-of-mass vector G, calculated in this section, stays unchanged, as Eq. (|7.14|) does not contribute to it. The 
Poincare algebra would not be fulfilled any more, therefore we have to drop the term proportional to Ja;' in Sec. IVl 
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APPENDIX A: POISSON-BRACKET ALGEBRA OF SPINNING PARTICLES STRESS-ENERGY 

TENSOR IN MINKOWSKI SPACE 

The following equal-time algebra, i.e., a;° = must be valid in Minkowski space [53| (see also [57|), 

l^mattcr(^)^ ^matter ^ _ [Wf^^^Xx) + Hr^'^'i^')] <5xx',, + d^^d.d'pd'g [/™„p,(x) (5xx'] , (Al) 
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l^matter(^)^ ^matter ^ _ 5xx',. " T„ (x') 5xx' j + d^O'^d'^ [g^npg{x) (5xx'] , (A2) 

l^mattcr(^)^ ^matter ^ _ ^^--""(.t) 5,,,,, - Hf^^^X^') ^xx' + d^d'g [h„,jq{x) 6^^,] . (A3) 

Here Jxx' = x'), where x and x' are the spatial parts of a; and x' . This local algebra is a consequence of the global 
Poincare algebra, whose generators can be written in terms of integrals over certain components of the stress-energy 
tensor, similar to Eq. p.8p . The terms containing /, g and h turn into vanishing surface terms in the generators of 
the Poincare algebra. It holds: 

fmnpq = fnmpq = fmnqpj the Same for g and h , (A4) 

fmnpq — fpqmn ; 

the same for h . (A5) 
An explicit calculation with the Minkowski versions of (|4.25p . (|4.23p and (|4.26p shows that we have to set 

fmnpq(x) = = ginpq(x) , (A6) 



hinjq (*^) 



Sq)(n'Pi)(j - S 



(np) (to — np) (np) {m — np) 



S, with P., ^5,,-^. (A7) 



Now the local algebra is fulfilled linear in the spin variables, as it should be, because our variables are known to be 
canonical in the Minkowski case. It was already noted in [67|, in the context of quantum field theory, that hinjq{x) 
does not generally vanish if fields with spin are present, in particular, spin-i fields. For fields with spin | one even 
has fmnpq 7^ 0, sce (6^. The consequences of nonvanishing hinjq{x) are considered in Appendix B. 



APPENDIX B: POISSON-BRACKET ALGEBRA OF NON-SPINNING PARTICLES STRESS-ENERGY 

TENSOR IN GENERAL RELATIVITY 

We assume that the following equal-time constraint algebra on the nonreduccd phase space without gauge fixing is 
valid HI,!!!!!!: 

{n{x),n{x')} = - [n^{x)Y'{x) + n^{x')Y'{x')] 6^^,,, , (bi) 
{n^ix),n{x')} = -nix)s^^,,^, (b2) 

{H^{x),Hjix')} = - Hj{x)5^^,^, - m{x') 5xx'j . (B3) 

Note that, compared to the algebra of the last section, the T^j term and the surface terms are absent. If this local 
algebra is fulfilled, then the global Poincare algebra can also be derived [13, 13 ■ At this point the constraints are not 
solved and no coordinate conditions are imposed, i.e., one has to use {7ij(x, t), 7r'^'(x', <)} — 167r(5fj(5(x — x'), where 

Remember that H and Hi are a sum of matter and field parts. The field-field Poisson-brackets cancel with the 
field terms on the right hand side of each relation of the algebra, because the algebra is fulfilled if no matter would 
be present, see [11]. In the context of (|2.4[) . the matter parts do not depend on tt*-', Eq. (|3.4[) . Therefore all terms 
proportional to tt*-' that arise from the mixed matter-field Poisson-brackets must vanish separately. These are exactly 
the ones with Ti^"^'^: 

{H^°^^(x),H™^"'='^(a;')} + {W'"''""(a;),7^*^'='^(a;')} = 0, (B4) 

These conditions are indeed equivalent to Eq. p.Sp . From (jBSP follows: 

^'='(x') f7fcKx07""(xO '^y"r^^^ - 2^^S?!^1 = (B6) 



Because of the factor 2, we have ^lyfci(x') = as the only possible solution. In Eq. (IB4|1 we make a general ansatz 



matter 



(57'J(x' 



^ = a.,(x)<5(x - x') + ^ a,/-'="(x)afc, . . . 9fc„5(x - x') , (B7) 



" 5(a.,...a.„7-) ' 
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integrate over x' and demand that the term with the highest number of derivatives on tt*-' must vanish separately. 
The resuhing equation is similar to (|B6p . this time leading to a^^'^'-'"'^" (x) = 0. Now we have effectively reduced N 

by one, proceeding this way we get a^j'^^'"''"{x.) — for all n with 1 < n < A^. Comparing with p.3p we see that aij 
has to be identified as ■^y/jTij. 

From (|3.5p then immediately follows: 

{nf^"'{x),n'^^''''^{x')} = ^T,u{x') [Sli'^'ix') S^^,j + j%{x') S^^] , (B9) 
{Hf ^(x),HJ"^""(a;')} = 0. (BIO) 
For the matter part therefore a similar algebra as in the Minkowski case has to hold, 

l^mattcr(^)^ -^matter (^/) I ^ _ [jif'''''''' {x)-f'^ (x) + Hf''''''' (x'h'^ (x')] (5xx'j , (Bll) 

{Hr""(a^),W""""(a;0} =-^"""'='■(a;)'5xx^.- V7T,fe(a;') [sh"' {x') S^^, .i + j%{x') 6^^,] , (B12) 
Hf^'^^^ix')} = - Hf^'^^^ix) 5xx'.» - Hf'''**" {x') S^^,,j . (B13) 

If the consistency conditions (|3.4p and (|3.5p hold, still in the context of (|2.4p . then this algebra can be used to 
validate the canonical variables of the matter part on the nonreduced phase space. This algebra is indeed fulfilled for 
point-masses, and of course also Eqs. (|2.4p . (|3.4p and (|3.5p . 

Because in the algebra (|B11||BT21) there arc no variations with respect to 7^ and tt'^ left any more, we can now 
consider its Minkowski space limit. This gives the algebra of the last Section with — fmnpqix) = ginpq{x) = hinjq{x). 
But in the last Section we have seen that for spinning bodies hinjq{x) does not vanish, not even in the Newtonian case. 
Therefore we already see by inspecting the Minkowski case that the coupling to gravity can not be of the simple kind 
defined by (|2.4p . p.4p and (|3.5p . Another problem that must be addressed by a gauge independent formulation is that 
we have (16. 2p and (|6.3p instead of Eq. p.Sp . The total divergence in (|6.2p contributes to the leading order. Therefore 
p.Sp is not even fulfilled at the leading order. This together with the nonvanishing hinjq{x) leads to additional 
contributions to the algebra (jBHIB3p and suggests its extension when spinning objects are present. Extensions to 
p.4p may also be considered, recall Ref. 55]. 

Important about the algebra of (first-class) constraints is their connection to the gauge structure of the considered 
theory, see, e.g., [6^. This makes the algebra (jBlHB3p quite robust even if other systems are coupled to gravity. 
Extended forms of the algebra (IB1HB3I) were, to the best of our knowledge, indeed only found when the gauge 
structure was also extended ^23. , 28„ 70. , ,71.] . This we will keep in mind when investigating higher approximations in 
future. 

The considerations of this Appendix do not show up any inconsistencies of the canonical formulation for spinning 
bodies given in this paper because they are requiring that no coordinate conditions are imposed. Rather they indicate 
that the gauge conditions (|2.9p must be seen as essential part of our formulation. 



APPENDIX C: PARTIE FINIE AND ANALYTIC REGULARIZATION 

In this Appendix we will present the regularization techniques we used in this work. We first give a short overview 
of Hadamard's partie finie method. Let us consider / being a real function defined in an environment of the point 
xq G M^, except in this point where / is singular. We define a family of complex- valued functions /„ as follows: 

U ■■ C3e ^ /„(£) = /(xo+en) eC. (CI) 

We expand /n in a Laurent-series around e = 0: 

oo 
m=-N 

The regularized value of / at Xq is defined as the coefficient at e*^ in the expansion (|C2p mean-valued over all unit 
vectors n, [I,!!!!!!, 

/,.eg(xo) = ^^dr!ao(n). (C3) 
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This formula can be used to calculate integrals with delta-distributions. We define 

d^xf{^)6{^-Za):^froe{Za), (C4) 

which provides us with a formula for calculating Poisson integrals of the form 



(C5) 



A complicated example is given by (|8.12p . 

Integrals that do not contain a delta function are regularized analytically [gT'I . First we perform all differentiations 
in the integrand, and then constrain ourselves to the two particle case. The integrand then depends on ri — |x — Zi|, 
rii = (x— Zi)/ri and r2, n2. Now we introduce the analytic regularization parameter e by replacing r" by r"'^'^'^ , and 
r2 by Tj^'^'^. The vectors n\ and are then rewritten as partial derivatives dj and dj with respect to the particle 
positions 

Qa a+l r . Q gaga a+2 

« = r»i---^ + ^T^, (C6) 

a + l a a[a + 2) 

(a-l)(a + l) (a-l)(a + l)(a + 3) ■ ^ ' 

These equations are sometimes not defined without regularization, because the right-hand side might diverge in special 
cases. The derivatives and dj ; 
[7^ to carry out the integrations: 



cases. The derivatives and are now pulled out in front of the integral, and we can use the famous formula from 



7 3 3./mr(^)r( 
J 'L. r(-f)r(-|)r(^ 



n^l^ ' ^/ { \ i^r"+'3+3. (C8) 



After the partial derivatives d} and dj with respect to the particle positions, that were pulled out of the integral 
before, are performed, the limit e ^ can be considered. In all cases emerging in our calculations this limit was 
independent of the direction in the (/i, j/)-plane. Astonishingly, the most complicated integral appearing in this work 
has the simplest solution: 

d3x/lJJ^f"0(2),,(/.(2),, =0. (C9) 



[1] R. Arnowitt, S. Deser, a nd C. W. M isner, in Gravitation: An Introduction to Current Research, edited by L. Witten (Wiley, 
New York 1962), p. 227; |arXiv:gr ^qc/0405109 

[2] T. Damour, P. Jaranowski, and G. Schafer, Phys. Rev. D 77, 064032 (2008). 

[3] J. Steinhoff, S. Hergt, and G. Schafer, Phys. Rev. D 77, 081501(R) (2008); arXiv:0712.1716. yl [gr-qc]. 

[4] T. Kimura, Prog. Theor. Phys. 26, 157 (1961). 

[5] T. Kimura and T. Toiya, Prog. Theor. Phys. 48, 316 (1972). 

[6] T. Ohta, H. Okamura, T. Kimura, and K. Hiida, Prog. Theor. Phys. 51, 1598 (1974). 

[7] T. Ohta, H. Okamura, and T. Kimura, Nuovo Cimento 27B, 103 (1975). 

[8] T. Ohta and T. Kimura, Prog. Theor. Phys. 81, 679 (1989). 

[9] G. Schafer, Ann. Phys. (N.Y.) 161, 81 (1985). 
[10] T. Damour and G. Schafer, Gen. Relativ. Gravit. 17, 879 (1985). 
[11] T. Damour and G. Schafer, Nuovo Cimento B 101, 127 (1988). 
[12] P. Jaranowski and G. Schafer, Phys. Rev. D 55, 4712 (1997). 
[13] P. Jaranowski and G. Schafer, Phys. Rev. D 57, 7274 (1998). 
[14] P. Jaranowski and G. Schafer, Phys. Rev. D 60, 124003 (1999). 
[15] T. Damour, P. Jaranowski, and G. Schafer, Phys. Lett. B 513, 147 (2001). 
[16] H. P. Kiinzle and J. M. Nester, J. Math. Phys. 25, 1009 (1984). 
[17] D. D. Holm, Physica 17D, 1 (1985). 



18 



[18] D. Bao, J. Marsden, and R. Walton, Comm. math. Phys. 99, 319 (1985). 

[19] L. Blanchet, T. Damour, and G. Schiifer, Mon. Not. R. Astr. Soc. 242, 289 (1990). 

[20] G. Schiifer, Astron. Nachr. 311, 213 (1990). 

[21] J. Schwinger, Phys. Rev. 130, 1253 (1963). 

[22] J. Schwinger, Phys. Rev. 132, 1317 (1963). 

[23] C. Teitelboim, in General Relativity and Gravitation. One Hundred Years After the Birth of Albert Einstein, Vol. 1, ed. A. 

Held (Plenum Press, New York 1980), p. 195. 
[24] P. A. M. Dirac, in Recent Developments in General Relativity (Pergamon Press, Oxford 1962), p. 191. 
[25] T. W. B. Kibble, J. Math. Phys. 4, 1433 (1963). 
[26] S. Deser and C. J. Isham, Phys. Rev. D 14, 2505 (1976). 
[27] J. Geheniau and M. Henneaux, Gen. Relativ. Gravit. 8, 611 (1977). 
[28] J. E. Nelson and C. Teitelboim, Ann. Phys. (N.Y.) 116, 86 (1978). 

[29] J. Isenberg and J. Nester, in General Relativity and Gravitation. One Hundred Years After the Birth of Albert Einstein, 

Vol. 1, ed. A. Held (Plenum Press, New York 1980), p. 23. 
[30] W. Tulczyjew, Acta Phys. Polon. 18, 37 and 393 (1959). 
[31] A. Papapetrou, Proc. Roy. Soc. Loud. 209, 248 (1951). 

[32] W. G. Dixon, in Isolated Gravitating Systems in General Relativity, Proceedings of the International School of Physics 

'Enrico Fermi', Varenna, Lake Como, ed. J. Ehlers (North-Holland, Amsterdam 1979), p. 156. 
[33] T. D. Newton and E. P. Wigner, Rev. Mod. Phys. 21, 400 (1949). 

[34] C. M0ller, Ann. Inst. H. Poincare 11, 251 (1949); The Theory of Relativity (Oxford Univ. Press, Oxford 1957). 

[35] E. Corinaldesi and A. Papapetrou, Proc. Roy. Soc. Lond. 209, 259 (1951). 

[36] A. D. Fokker, Relativiteitstheorie (P. Noordhoff, Groningen 1929). 

[37] M. H. L. Pryce, Proc. Roy. Soc. Lond. 195A, 62 (1948). 

[38] J. L. Synge, Phys. Rev. 47, 760 (1935). 

[39] G. N. Fleming, Phys. Rev. 137, B 188 (1965). 

[40] C. Lanczos, Z. Phys. 59, 514 (1929). 

[41] M. Mathisson, Acta Phys. Polon. 6, 163 (1937). 

[42] F. A. E. Pirani, Acta Phys. Polon. 15, 389 (1956). 

[43] O. Semerak, Mon. Not. R. Astron. Soc. 308, 863 (1999). 

[44] L. E. Kidder, C. M. Will, and A. G. Wiseman, Phys. Rev. D 47, R4183 (1993); L. E. Kidder, Phys. Rev. D 52, 821 (1995). 
[45] H. Tagoshi, A. Ohashi, and B. J. Owen, Phys. Rev. D 63, 044006 (2001); G. Faye, L. Blanchet, and A. Buonanno, Phys. 

Rev. D 74, 104033 (2006); L. Blanchet, A. Buonanno, and G. Faye, Phys. Rev. D 74, 104034 (2006); 75, 049903(E) (2007). 
[46] W. D. Goldberger and I. Z. Rothstein, Gen. Relativ. Gravit. 38, 1537 (2006); R. A. Porto, Phys. Rev. D 73, 104031 (2006). 
[47] T. Regge and C. Teitelboim, Ann. Phys. (N.Y.) 88, 286 (1974). 

[48] A. Hanson, T. Regge, and C. Teitelboim, Constrained Hamiltonian Systems (Academia Nazionale dei Lincei, Roma 1976). 
[49] S. Hergt and G. Schiifer, Phys. Rev. D 77, 104001 (2008). 
[50] K. Yee and M. Bander, Phys. Rev. D 48, 2797 (1993). 
[51] H. P. Kiinzle, J. Math. Phys. 13, 739 (1972). 

[52] J. Natario, Commun. Math. Phys. 281, 387 (2008); [axXiv:gr-qc/070"308T| 

[53] B. S. DeWitt, Phys. Rev. 160, 1113 (1967). 

[54] D. G. Boulware and S. Deser, J. Math. Phys. 8, 1468 (1967). 

[55] K. Kuchaf, J. Math. Phys. 18, 1589 (1977). 

[56] A. Trautman, Gen. Relativ. Gravit. 34, 721 (2002). 

[57] J. Schwinger, Particles, Sources, and Fields, Vol. I (Perseus Books, Reading/Massachusetts 1998). 

[58] A. Ohashi, Phys. Rev. D 68, 044009 (2003). 

[59] J. M. Martm-Garci'a, see |http://metric.iem.csic.es/Martin-Garcia/xAct/index.html[ 

[60] S. Wolfram, Mathematica. A System for Doing Mathematics by Computer, 2nd ed. (Addison- Wesley, Redwood City 1991). 
[61] P. Jaranowski, in Mathematics of Gravitation, Part H: Gravitational Wave Detection, ed. A. Krolak (Banach Center 

Publications, Vol. 41, Part II, Warszawa 1997), p. 55. 
[62] L. Blanchet and G. Faye, J. Math. Phys. 41, 7675 (2000). 
[63] R. A. Porto and I. Z. Rothstein, Phys. Rev. Lett. 97, 021101 (2006). 

[64] R. A. Porto and I. Z. Rothstein, Comment on 'On the next-to-leading order gravitational spin(l)-spin(2) dynamics' by J. 

Steinhoff et al. larXiv:0712.2032V l [gr-qc]. 
[65] R. A. Porto, New results at 3PN via an effective field theory of gravity, Proc. of the 11th Marcel Grossmann Meeting, 

Berlin, July 23 -25, 2006; arXiv;gr-qc/0701106vl. 
[66] T. Damour, P. Jaranowski, and G. Schiifer, Phys. Rev. D 62, 021501(R) (2000); 63, 029903(E) (2001). 
[67] J. Schwinger, Phys. Rev. 127, 324 (1962). 
[68] J. Schwinger, Phys. Rev. 130, 800 (1963). 
[69] L. Castellani, Ann. Phys. (N.Y.) 143, 357 (1982). 
[70] M. Henneaux, Phys. Rev. D 27, 986 (1983). 

[71] L. Castellani, P. van Nieuwenhuizen, and M. Pilati, Phys. Rev. D 26, 352 (1982). 
[72] M. Riesz, Acta Math. 81, 1 (1949). 



